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ABSTRACT 

The observation by the Swift X-ray Telescope of the Fe Kai,a2 doublet during a large flare on 
the RS CVn binary system II Peg represents one of only two firm detections to date of photosphcric 
Fe Ka from a star other than our Sun. We present models of the Fe Ka equivalent widths reported 
in the literature for the II Peg observations and show that they are most probably due to fluorescence 
following inner shell photoionisation of quasi-neutral Fe by the flare X-rays. Our models constrain 
the maximum height of flare the to 0.15 R» assuming solar abundances for the photospheric material, 
and 0.1 R* and 0.06 R* assuming depleted photospheric abundances ([M/H] = -0.2 and [M/H] = -0.4, 
respectively). Accounting for an extended loop geometry has the effect of increasing the estimated 
flare heights by a factor of ~3. These predictions are consistent with those derived using results of 
flaring loop models, which are also used to estimate the flaring loop properties and energetics. From 
loop models we estimate a flare loop height of 0.13 R„, plasma density of ~ 4 x 10 12 cm -3 and emitting 
volume of ~ 6 x 10 30 cm 3 . Our estimates for the flare dimensions and density allow us to estimate the 
conductive energy losses to E con d < 2 x 10 36 erg, consistent with upper limits previously obtained in 
the literature. Finally, we estimate the average energy output of this large flare to be ~ 10 33 erg sec -1 , 
or l/10th of the stellar bolometric luminosity. 
Subject headings: X-rays: stars — stars: coronae — stars: individual (II Pegasi) — stars: flare 



1. INTRODUCTION 

A fraction of the X-rays emitted during stellar flares 
or by hot stellar coronae are directed downwards towards 
the stellar photosphere. Here they interact with the pho- 
tospheric gas and are reprocessed through scattering and 
photoionisation events. The downward cascade following 
inner shell photoionisation of quasi-neutral gas in the 
stellar photosphere produces characteristic fluorescent 
emission from astrophysically abundant species. While 
the reprocessed X-ray spectra reach the observer at low 
flux levels compared to that of the flares/coronae, fluo- 
rescent emission lines, can be significantly stronger and 
can be detected against the spectra of the flare/corona 
itself. 

The intensity or equivalent width of the fluorescent 
lines detected at Earth depends uniquely on the details 
of the fluorescing spectrum, the photospheric abundance 
of the photoionised species and the geometry of the sys- 
tem (flare height and inclination angle; e.g. see fig. 2 
of Testa et al., 2008), making them potentially power- 
ful diagnostics. A number of theoretical studies exist 
which provide a framework for the interpretation of the 
strong iron fluorescent emission in the solar context (e.g. 
Tomblin, 1972; Bai, 1979), and more recently for ar- 
bitrarily photoionised slabs (Kallman et al., 2004) and 
stellar photospheres (Drake, Ercolano and Swartz, 2008, 
DES08). Drake & Ercolano (2008a,b) have also recently 
showed that fluorescent emission from neon and oxygen 



should also be detectable in the solar spectrum with cur- 
rent and future instrumentation. These may provide an 
independent measure of the photospheric solar oxygen 
and neon abundance, elements that are central to the so- 
lution of the "solar oxygen crisis" (e.g. Ayres et al. 2006, 
Socas-Navarro 2007, Basu & Antia 2008). 

Observationally, the 2s — \p 6.4 KeV Fe Ka doublet 
from low ionisation stages of iron has often been detected 
in solar spectra (e.g. Neupert et al., 1967; Doschek et al., 
1971; Fedelman et al., 1980; Tanaka et al., 1984; Parmar 
et al. 1984; Zarro et al., 1992), from disks around pre- 
main sequence stars (e.g. Tsujimoto et al., 2005; Favata 
et al., 2005), and from the photospheres of the single G- 
type giant HR 9024 (Testa et al., 2007) and of the RS 
CVn binary system II Peg (Osten et al., 2007, hereafter 
O07). 

The detection of the 6.4 keV fluorescent iron line in the 
Swift X-ray Telescope spectrum of II Peg taken during 
a large flare (O07) is particularly interesting as it rep- 
resents the first of only two detections (with HR 9024, 
Testa et al., 2007) of photospheric fluorescence emission 
in stars other than the Sun. O07, however, assigned the 
excitation mechanism to electron impact ionization of 
photospheric Fe, rather than photoionisation. In this pa- 
per we show by means of Monte Carlo fluorescence mod- 
eling that the Fe Ka equivalent widths reported by O07 
are perfectly consistent with photoionisation induced flu- 
orescence. We also show that our predictions of the flar- 
ing loop scale height agree with those from flaring loop 
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Fig. 1. — Example Fc Ka equivalent widths as a function of flare height and inclination angle for three different photospheric metallicities 
for the la segment of the flare, as defined by O07. The observed value and its error is indicated by the horizontal black dashed line and 
teh shaded area. Left: Models calculated at solar abundances (Grevesse & Sauval, 1998); Central: Models calculated with [M/H] = -0.2; 
Right: Models calculated with [M/H] = -0.4. 



models (Reale, 2007) and discuss other implications on 
flaring plasma density, volume and energetics. 

In Section 2 we briefly describe our Monte Carlo model, 
and in Section 3 we present and discuss our results. Sec- 
tion 4 deals with the flaring loop properties and energet- 
ics from hydrodynamic models. Our final conclusions are 
given in Section 5. 

2. MONTE CARLO FLUORESCENCE CALCULATIONS 

The 3D Monte Carlo calculations presented here were 
performed using a modified version of the photoionisation 
and dust radiative transfer code, MOCASSIN (Ercolano 
et al. 2003, 2005, 2008). The code uses a a stochas- 
tic approach to the radiation transfer allowing it to deal 
with problems of arbitrary geometry, while treating both 
the primary and secondary components of the radiation 
field self-consistently. The modified version used in this 
work, which can deal with the production and transfer 
of fluorescent radiation, is further described by DES08. 

We performed calculations to fit all Fe Ka detections 
reported by O07 for the 5 different orbits (flare inter- 
vals) of the Swift X-ray telescope. Spectra representing 
the flaring source of X-ray irradiation for input into the 
photospheric fluorescence modeling were computed for 
each of these flare intervals. Spectra were computed us- 
ing the CHIANTI database version 5.2 (Dere et al. 1997, 
Landi et al. 2006) for the energy range 6.2-30 keV, using 
the temperatures and emission measures found by O07 
from two-temperature model fitting and listed in Table 3 
of that article. The upper limit to the energy range was 
chosen to be sufficiently high that any contribution to the 
observed fluorescent Fe K line from higher energy pho- 
tons in these thermal spectra should be less than 1%. 

Photospheric metal abundances in II Peg have been 
studied by Ottmann et al. (1998) and Berdyugina et al. 
(1998). Both are based on high resolution, high signal-to- 
noise spectroscopy, and there is little to choose between 
the analyses. The former obtained abundances [Fe/H] = 
-0.2, [Mg/H] = -0.15, and [Si/Hj = -0.15, expressed 
relative to the solar composition in conventional spec- 
troscopic logarithmic bracket notation, and each with 
uncertainty of ±0.1. As a verification of their method, 
they successfully recovered the currently accepted solar 
atmospheric parameters from a spectrum of moonlight. 
Berdyugina et al. (1998) obtained [M/H] = -0.4 ± 0.1, 
and verified their approach against the detailed analy- 
sis of the K0 giant Pollux by Drake & Smith (1991). 
It is therefore reasonable to infer that II Peg is mildly 
metal-poor by 0.1-0.4 dex relative to the Sun. We have 
investigated the effects of this slight metal depletion on 
Fe Ka equivalent widths by running three sets of models: 



(i) with solar abundances (Grevesse & Sauval, 1998), (ii) 
with [M/H] = -0.2 and (iii) with [M/H]= -0.4. 

Fe Ka equivalent widths also strongly depend on the 
geometry of the system, i.e. height of the flare, h, and 
inclination angle, 9. In particular Fe Ka efficiencies 
(and therefore equivalent widths) depend on the function 
f(h,0), which has been numerically determined for the 
case of a single flare illuminating a spherical photosphere 
by DES08. Assuming that the flare can be approximated 
by a point source above the stellar surface, we have per- 
formed calculations at ft, = 0. ,0.1, 0.2, 0.3, 0.5 and 9 = 0. 
and used the analytical form of f(9) given by DES08 in 
order to estimate the inclination angle dependence. We 
note that in order to scale 9 = equivalent widths to an 
arbitrary value of 9 using the f(6) functions one has to 
first scale the latter such that /(0) = 1. 

3. RESULTS FROM MONTE CARLO FLUORESCENCE 
MODELING 

3.1. Assuming a point source flare 

Representative flare geometry diagnostic diagrams il- 
lustrating the Fe Ka line equivalent width as a function 
of 9 for different flare heights are presented in Figure 1 for 
the three values of [M/H] (0.0, -0.2 and -0.4) for Orbit 
segment la. Also indicated in these figures for compar- 
ison is the observed equivalent width (dashed line) and 
uncertainties (shaded area). The combination of pre- 
dicted and observed equivalent widths can be combined 
so as to determine the inferred flare height implied by the 
observations as a function of the angle 9. For the grid of 
calculations comprising the five different flare segments 
and three values of metallicity, we have recast the results 
into this form in which the flare point source height im- 
plied by the observed equivalent width is illustrated as a 
function of 9 in Figure 2. 

At solar metallicities, Fe Ka equivalent widths re- 
ported for Orbits la, lb, lc and 2b are well-fit by mod- 
els with flare height, h, up to ^0.15 R* for low (<20°) 
inclination angles; larger values of 9 would necessarily 
demand a lower h. Lower metallicity models can also 
successfully reproduce the observations of Orbits la, lb, 
lc and 2b. Models with [M/H]= -0.2 (central column) 
indicate a maximum flare height of ~ 0.1 R*, while mod- 
els with [M/H]= —0.4 indicate a maximum flare height of 
~0.06 R* . We note that the lowest metallicity models fall 
slightly short of reproducing the equivalent width of the 
Fe Ka line measured during Orbit 2b. The discrepancy 
is small enough to be easily explained by observational 
uncertainties in the measured equivalent width, which 
appears somewhat larger compared with the values de- 
termined for the other time intervals by O07. 
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Finally, the observations during Orbit 2a yield much 
lower values of Fe Ka equivalent widths (by a factor of 
~3), than for Orbits la, lb, lc and 2b. This was already 
pointed out by O07, although there was no obvious ex- 
planation for the difference. According to our model the 
Fe Ka equivalent width observed during Orbit 2a implies 
maximum flare heights between 0.5 (for [M/H] = 0) and 
0.45 R* (for [M/H] = -0.4) with an uncertainty of a factor 
of ~ 2. These values are even slightly discrepant with the 
those deduced from the other flare intervals. Intrigued 
by this discrepancy, we have independently verified from 
the Swift data that the line is indeed weak during this 
segment, though we obtain a slightly larger 95% confi- 
dence upper limit of EW < 43 eV than found by O07. 

3.2. Effects of extended loop geometry on Fe Ka 
fluorescence 

The above fluorescence calculations and interpretation 
assume that the flare can be approximated by a point 
source of X-rays. In the case of large flares similar to 
those seen on the Sun, with peak temperatures of order 
10 7 K, the "hard" X-ray photons above the 7.11 keV Fe K 
ionisation threshold will be emitted almost uniquely by 
the hotter loop apex and the point source approxima- 
tion should be accurate. However, Testa et al. (2008) 
have shown that the extreme flare observed on the giant 
HR 9024 that reached temperatures of ~ 10 8 K emits 
ionising radiation over its entire length, such that the 
equivalent point-source flare height for Fe Ka fluores- 
cence was /i/3, where h is the actual loop height. 

The II Peg flare analysed here is similar to the HR 9024 
one in terms of plasma temperature, and we can expect 
a similar scaling of the point-source height and true flare 
height. The corresponding heights for the [M/H] = —0.2 
and —0.4 models are h ~ 0.3 and 0.15 R*, respectively, 
taking segments la, lb and lc as representative. 

4. FLARING LOOP PROPERTIES AND ENERGETICS 

The allowed ranges of flaring loop scale height derived 
from the strength of the Fe K line in the previous section 
can be compared with predictions from detailed hydro- 
dynamic flare models. Here, we use the scaling laws for 
loop length, cross-sectional area and plasma density de- 
rived from hydrodynamic simulations by Reale (2007) in 
application to the light curve and plasma emission mea- 
sure and temperature found from the Swift observations 
by O07. 

Relevant quantities required by the Reale (2007) rela- 
tions are the flare rise time ijv/,3 (in units of 10 3 s), the 
maximum plasma temperature Tq^ (in units of 10 7 K), 
the peak emission measure $ (cm -3 ), and the 1/e decay 
time t' s (s). We estimate the flare rise and decay times 
from Fig. 1 in O07 al. (2007), and take the flare tem- 
perature and emission measure data from their Table 2. 
Since we are interested in the peak flare temperature, 
we adopt the highest temperature obtained from their 
3-temperature best-fit models. The final set of adopted 
parameters is 
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Fig. 2. — Inferred flare height vs astrocentric angle for the differ- 
ent flare segments computed for the different photospheric metal- 
licities considered here: [M/H]=0.0 (top); [M/H]= -0.2 (middle); 
and [M/H]= —0.4 (bottom). The observed X-ray spectral proper- 
ties and Fe Ka equivalent widths for segments lb and lc are es- 
sentially identical and have been averaged. The different segments 
are colour-coded, with upper and lower limits given by dashed and 
dotted curves, respectively. 

Here, we have assumed that the maximum plasma tem- 
perature is the value observed at the time of maximum 
density, and that this time corresponds to the peak in 
the observed flare emission. Since the flaring loop is a 
closed volume, the latter assumption is straightforward. 
However, as noted by Reale (2007), the maximum flare 
temperature can be reached some time after the peak 
density is reached; we allow for this in the analysis be- 
low. 

4.1. Loop length 

The flaring loop half-length can be estimated from the 
from rise time and the maximum plasma temperature, 
as given in Eqn. (12) of Reale (2007): 



U « 3 ^ 2 T o yt M ,3 (5) 

where Lg is the loop half-length in units of 10 9 cm. The 
parameter ip is the ratio of the maximum plasma tem- 
perature to that at density maximum, 
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From Fig. 2 of O07, it can be seen that the rise phase 
flattens toward the emission maximum; this suggests the 
flaring loop is close to equilibrium at the maximum and 
that, therefore, ip ss 1 and the assumption of maximum 
plasma temperature at the time of maximum density 
holds. To allow for deviations from this assumption, we 
also quote results for -0 « 1.3 (see examples given in 
Reale 2007). For T , 7 w 30, t M ,3 «2>«1 

L 9 w30 
For r/> rs 1.3 

L 9 «70 
As a reasonable average value, we adopt Lg s=a 50 ± 20. 
Adopting a radius for II Peg of 3.4i? Q (Berdyugina et al. 
1998), this loop half-length corresponds to a loop height 
of ~ 0.13i?*. This compares very favourably with the 
values we have estimated from the Fe K fluorescence line 
in the previous section, provided the flare were located 
toward the centre of the visible stellar disk and not to- 
ward the limb. 

As a verification that this loop length value is reason- 
able, we can also obtain an estimate of the upper limit 
to Lg from the decay timescale. If the decay is entirely 
driven by cooling, with no further heat input, the ob- 
served decay timescale can be equated with the thermo- 
dynamic decay time, 

r' s = 4>t s (6) 

where <6 « 1.3 and 



3.7 x 10~ 4 



= 120- 



U 



(7) 



For t' s w 9000 s, we obtain an upper limit for the loop 
length: 

L g < 320. 

While much larger, this upper limit is consistent with the 
value obtained from the rise phase above. 

Hucnemocrdcr et al. (2001) observed a large flare on 
II Peg during a 45 ks observation with the Chandra High 
Energy Transmission Grating Spectrometer (HETGS), 
during which the count rate was observed to rise to 2.5 
times its quiescent value and the peak flare temperature 
reached logT ^7.6. While much less intense than the 
Swift flare discussed here, the decay timescale of ~ 65 ks 
was somewhat longer. Hucndmoerder et al. (2001) in- 
terpreted the light curve in terms of both single loop 
and two-ribbon type flares and inferred a loop height 
in the range 0.05-0.25 R* for the former, depending on 
the plasma density assumed. This height then scales in- 
versely with the cube of the number of loops involved — 
e.g. smaller by a factor of 5 for 100 loop strands. While 
this estimate must be considered very approximate, the 
similarity with the loop height derived here for the Swift 
flare suggests that both events occurred in similar coronal 
structures that might be typical of the outer atmosphere 
of II Peg. 

4.2. Plasma Density, Volume and Loop Cross- sectional 

Area 

From Reale (2007), the maximum possible density 
equilibrium value at the loop apex in units of 10 10 cm~ 3 
is given by 



«o,io = 13 



J 0,7 
La 



While this is an upper limit, the true value should not 
be lower by more than a factor of 2. For Toj s=s 30 and 
Lg rs 50, the density is then 

^o,io < 360 
A reasonable value to adopt is no,io ~ 200. This value 
is compatible with density diagnostics from line ratios 
of hot He-like triplets, such as Mg XI, observed during 
flaring and quiescence on II Peg (Hucncmoerder et al. 
2001; Testa, Drake & Peres 2004). Similarly high den- 
sities have also been reported by Gudel et al. (2002) 
during an X-ray flare in Proxima Centauri. From Reale 
(2007), the volume is: 

v^-r- (8) 



where 
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Taking T avg rs 150 MK from the single temperature 
model fit of O07 (their Table 2), we obtain 



4x 10 1 



and 



U « 6 x 10 3 



The cross-sectional area is: 
A 



V 



(10) 



For Lg « 50, A rj 6 x 10 19 cm 2 . For a single loop with cir- 
cular cross-section we obtain a radius in units of 10 9 cm 
of rg « 4, which is slightly less than 1/10 the loop half- 
length and very similar to the values found for solar flar- 
ing loops (e.g. Cheng et al. 1980, Golub et al. 1980, Peres 
et al. 1987). 

4.3. Conductive energy losses 

One of the most striking aspects of the Swift II Peg 
flare is the potentially large flare energy budget. Us- 
ing the classical Spitzer formulism, O07 found an upper 
limit to the conductive losses of 5 x 10 43 erg, but were 
unable to constrain this further without estimates of the 
flare dimensions and density. The bolomctric luminos- 
ity of II Peg is log L ho i = 34.2 (Marino et al. 1999), and 
conductive losses as high as 10 43 erg during a single, com- 
pact flare would represent an astounding concentration 
of energy. 

We can refine the estimate of the conductive energy 
losses using Eqn. (5) in O07 and the parameters derived 
above for the flaring loop: 

_ nT 7 / 2 <i>At 

Econd - L 2 n 2 . 

where, for Orbit 1, T w 300 MK (an upper limit for the 
flare as a whole), $ (=s 10 56 , k rs 10~ 6 (in c.g.s. units), 

At » 2 ks, Lg rj 50, n e w n avg ~ ' : "< ! ~ ' 

obtain 



E n 



< 2 x 10 



:!(> 



erg 



which is similar to the total radiated energy of 6 x 10 36 erg 
in the 0.01-200 keV energy band estimated by O07. The 
duration of the flare was of order 10 4 s, and the radia- 
tive and conductive losses therefore represent an energy 
output of ~ 10 33 erg s _1 , or about 1/10 the stellar bolo- 
mctric luminosity. 
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5. CONCLUSIONS 

We have presented a set of Monte Carlo calculations 
able to reproduce the Fe Ka equivalent widths from the 
Swift X-ray Telescope observations of II Peg during a 
'superfiare' (Osten et al., 2007, O07). Our models show 
that the data are consistent with the Fe Ka emission 
being produced by fluorescence following K-shell pho- 
toionisation of quasi-neutral iron in the stellar photo- 
sphere. This contrasts with the interpretation of O07, 
who favoured collisional ionisation by non-thermal elec- 
trons to produce the Fe Ka emission. They argued that 
the normal incidence flare X-ray penetration depth re- 
quired in the photosphere to obtain the observed equiv- 
alent widths is similar to the r = 1 Compton scatter- 
ing depth, and consequently any fluorescent photons pro- 
duced would not easily escape. However, this assessment 
was based on a simple analytical formula appropriate for 
optically thin cases in which only a small fraction of the 
incident X-ray flux undergoes photoabsorption or scat- 
tering (e.g. Liedahl, 1999; Krolik & Kallman, 1987). The 
semi-infinite photospheric case lies outside the range of 
applicability of this formula. Furthermore, incident an- 
gles on the photosphere range from ~ 0-90° and path 
lengths for escape can therefore be much smaller than gas 
penetration depths (by a factor equal to the inverse co- 
sine of these angles) . Our stochastic treatment includes 
Compton scattering and shows that the observed equiv- 
alent widths can be produced by fluorescence. We also 
note that the impact excitation mechanism is a very low 
efficiency process and requires a large amount of energy 
in the form of accelerated electrons. This was already 
noted by, e.g., Parmar (1984) in reference to solar Fe flu- 
orescence and by Ballantyne (2003), and was discussed in 
reference to a large flare on II Peg by Testa et al. (2008). 

The derived flare loop heights are h = 0.15 R», assum- 
ing solar photospheric abundances, h = 0.1 R* assuming 
photospheric abundances depleted by [M/H] = -0.2 and 
h = 0.06 R* assuming [M/H] = -0.4. These value are in 
good agreement with the predictions from our alterna- 
tive analysis based on hydrodynamic models (e.g. Reale, 
2007) which yield h~ 0.13 R*. We estimate flaring loop 
properties and energetics using scaling laws based on the 
hydrodynamic loop modeling of Reale (2007) and ob- 
tain a plasma density of 4 x 10 12 cm~ 3 and volume of 
6 x 10 30 cm 3 . Using the derived values for the flare di- 
mensions and densities we obtain a more stringent upper 
limit for the conductive energy loss, E con d, than previ- 
ously possible, setting E con d < 2 x 10 36 erg. Considering 
the duration of the flare and the radiative and conduc- 
tive losses, the energy output of the flare is estimated at 
~10 33 erg sec -1 , which is approximately one tenth of the 
stellar bolometric luminosity. 

While we cannot rule out a contribution to the ob- 
served Fe Ka flux from electron impact with nonthermal 
electrons, as proposed by O07, such contribution is not 
necessary to explain the available data. 
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